A Black Hole in the Galactic Center Complex IRS 13E? 
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ABSTRACT 

The IRS 13E complex is an unusual concentration of massive, early-type stars 
at a projected distance of ~0.13pc from the Milky Way's central supermassive 
black hole Sagittarius A* (Sgr A*). Because of their similar proper motion and 
their common nature as massive, young stars it has recently been suggested 
that IRS 13E may be the remnant of a massive stellar cluster containing an 
intermediate-mass black hole (IMBH) that binds its members gravitationally in 
the tidal field of Sgr A*. Here, we present an analysis of the proper motions in the 
IRS 13E environment that combines the currently best available data with a time 
line of 10 years^. We find that an IMBH in IRS 13E must have a minimum mass 
of ~lO^M0in order to bind the source complex gravitationally. This high mass 
limit in combination with the absence so far of compelling evidence for a non- 
thermal radio and X-ray source in IRS 13E make it appear unlikely that an IMBH 
exists in IRS 13E that is sufficiently massive to bind the system gravitationally. 

Subject headings: Galaxy: center - Galaxy: nucleus - infrared: stars 
^Also: Department of Physics, University of California, Berkeley, CA 94720 

^Based on observations at the Very Large Telescope (VLT) of the European Southern Observatory (ESO) 
on Paranal in Chile 
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1. Introduction 

In spite of factors such as a strong tidal field due to the central supermassive black 
hole Sgr A* (e.g., Schodel et al. 2002; Ghez et al. 2003), or strong stellar winds, which 
pose serious obstacles for star formation, surprisingly, numerous young (a few times lO^yrs), 
massive stars can be found in the central parsec, e.g., the one to two dozen bright so-called 
He-stars (with characteristic He I emission lines, see Krabbe et al. 1995; Paumard et al. 
2001). The young, massive stars are mainly concentrated in the region ~10" in projection 
around Sgr A* and appear to rotate around Sgr A* in two coTinter-rotating disks that contain 
an apparently coeval population of Wolf-Rayet (WR), Luminous Blue Variables (LBV), and 
0/B-stars (Genzel et al. 2003; Levin & Beloborodov 2003, also Paumard et al. 2005, in 
preparation). Also, there are B-type main sequence stars in the immediate vicinity of Sgr A* 
(e.g., Ghez et al. 2003; Eisenhauer et al. 2005). The presence of the young, massive stars 
near Sgr A* is not well understood and various explanations for their existence are currently 
being discussed, such as, e.g., infall and coUision of molecular clouds, infall and dissolution 
of a massive cluster (e.g., see discussions in Genzel et al. 2003), or star formation in a 
self-gravitating accretion disk (e.g., Nayakshin & Cuadra 2004; Milosavljevic & Loeb 2004). 

In this context, the IRS 13E complex is a very intriguing object. Located ~ 3.5" in 
projection from Sgr A*, it comprises almost half a dozen massive stars within a projected 
radius of ~0.25". These stars appear to be WR stars or 0-type supergiants (e.g., Krabbe 
et al. 1995; Eckart et al. 2004; Maillard et al. 2004). In addition to representing an unusual 
concentration of bright, early-type stars, the IRS 13E sources also have very similar proper 
motions (see, e.g.. Fig. 18 in Genzel et al. 2003). Maillard et al. (2004) suggested that 
the IRS 13E complex is the remnant core of a massive star cluster that had fallen into the 
central parsec and dissolved there. Their main arguments are a) the similar proper motions 
of the IRS 13E cluster members and b) their finding that the stars in IRS 13E appear to 
be massive, short-lived, and therefore young. Since one would expect that an association 
of stars such as IRS 13E should be disrupted by the tidal field of Sgr A*, they speculated 
that an intermediate-mass black hole (IMBH) in IRS 13E may prevent its disruption. An 
IMBH in the core of an infalling cluster would also provide an effective means of inward 
transportation through increased dynamical friction as has been suggested by Hansen & 
Milosavljevic (2003) (see, however, Kim, Figer, & Morris 2004). 

The possible existence of an IMBH in the GC is currently a hypothesis of high interest 
in the field. Therefore we consider it timely to present the newest data on stellar dynamics 
in IRS 13E. In this letter, we analyze proper motions of the stars in and near the IRS 13E 
complex, derived from ten years of near-infrared (NIR) speckle and AO observations of the 
GC stellar cluster. 
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2. Observations and Data Reduction 

Observations of the GC stellar cluster are routinely performed since spring 2002 with the 
CONICA/NAOS NIR camera/adaptive optics (AO) system at the ESO VLT unit telescope 4 
on Cerro Paranal in Chile. For this work, we used K-band (and some H-band) imaging data. 
After standard data reduction (sky subtraction, dead pixel correction, flat fielding, shift- 
and-add) the final images were Lucy-Richardson (LR) deconvolved and beam restored with 
a Gaussian beam of ~60mas (40 mas for H-band images, respectively), corresponding to 
the diffraction limit of the VLT at 2.2 /xm. Stellar positions were extracted with StarFinder 
(Diolaiti et al. 2000). They were transformed into a coordinate system relative to Sgr A* 
with the aid of 18 reference stars with well known positions and proper motions, taken from 
(Ott 2004). Each of the data sets was divided into two parts that were reduced and analyzed 
independently. Uncertainties of the source positions were thus obtained by a comparison of 
two independent measurements for each epoch. 

We added to these data Gemini North AO H and K imaging data from July 2000. 
The images are part of the Gemini North Galactic Center Demonstration Science Data Set. 
Both images were LR deconvolved and beam restored. Source positions were extracted in 
the way described above (using the two independent measurements provided by the H and 
K images). In order to obtain a larger time baseline, we also used SHARP speckle imaging 
data from July 1995, June 1996, and June 2001. Details on the observation and reduction 
of SHARP imaging data may be found, e.g., in Eckart & Genzel (1997). Again, we applied 
LR deconvolution and used two independent data sets for each epoch. 

3. Dynamics in the IRS 13E complex 

The IRS 13E sources are named in Fig. 1. Eckart et al. (2004) label E3 and E4 as E3c 
and E3N, respectively. Maillard et al. (2004) noted that E3 is a multiple source and named 
the two components E3A and E3B. On the image presented here, E3 appears to have more 
than two components. We just label the brightest one of them E3 (the image in Fig. 1 is 
based on data with a Strehl > 40% in contrast to the Gemini image of Strehl ~5% used 
by Maillard et al. (2004)). As concerns E5, it appears highly confused with other sources 
and/or extended. 

Proper motions of the stars were determined by linear fits to the measured time- 
dependent positions with their respective uncertainties (adopting a GC distance of 7.9 kpc 
from Eisenhauer et al. 2003). In Fig. 1, the derived proper motions of all stars in and near 
IRS 13E are shown superposed on an image. The common proper motion of the IRS 13E 
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cluster members stands out clearly. We identify five stars within a region of about 0.5" 
diameter that could be part of IRS 13E (see also Table 1). There may be further poten- 
tial members, but they are either too weak, embedded in extended emission, or too close to 
brighter stars, such that no reliable proper motions could be determined. The proper motion 
of E6, while pointing into the same direction as the ones of the other IRS 13E sources, is 
significantly smaller in magnitude. Also, Maillard et al. (2004) note that E6 is extincted 
much less than the other IRS 13E stars. It is therefore questionable whether E6 forms part 
of IRS 13E and we omit it from the following analysis. The proper motion of E5 could not 
be reliably determined. 

Maillard et al. (2004) considered that IRS 13E might be bound by an IMBH. From spec- 
troscopically estimated radial velocities for E2 and E4 and assuming that the hypothetical 
IMBH is located halfway between the two stars, they estimated its mass to ~1000 Mq. With 
the proper motions, we can derive the minimum mass that is required to bind the IRS 13E 
cluster members. Here, it is first necessary to subtract the systemic motion of IRS 13E. 
Subsequently, the residual velocities of the individual sources can be used for estimating the 
mass in IRS 13E under the assumption that it is a gravitationally bound system. Since 
we cannot know the true systemic velocity, we examined three ) Using the average 

proper motion of the ensemble, b) using the proper motion of El, and c) using the proper 
motion of E4. Thus, all realistic cases, even extreme ones, should be covered. In the follow- 
ing, we limit our analysis to case a), which is the most conservative one because it results in 
the lowest enclosed masses. 

The average proper motion of IRS 13E, derived from the mean of the velocities of 
sources El, E2, E3, and E4, is 248 ± 25kms~^ westwards and 80 ± 48kms~-^ northwards. 
Here, we chose an unweighted average because we cannot know which one of the individual 
proper motions corresponds best to the possible systemic velocity. The proper motions of 
the IRS 13E sources after subtraction of a systemic velocity corresponding to their mean 
motion are shown in the left panel of Fig. 2. 

A first rough estimate of the mass needed to bind IRS 13E can be obtained from the 
residual velocities, through M — {v'^)R/G, where G is the gravitational constant and R 
lOmpc the size of the system. With (vj^j^) ~ 50kms^^((i'|)gp) ~ lOOkms"^) we obtain 
M Ki 5600 M0(M fa 22500 M0). 

We chose two more refined approaches to estimating the mass of a putative IMBH that 
gravitationally binds IRS 13E: a) Estimating the mass in IRS 13E with the Leonard-Merritt 
(LM) mass estimator (eq. 19 in Leonard & Merritt 1989) and b) deriving a lower hmit on 

the mass through the relation < = M, with R being the projected distance of a 

star from the IMBH, Vpj-oj its velocity projected onto the plane of the sky, r its real distance 
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from the IMBH, Vesc its escape velocity, G the gravitational constant, and M the mass of 
the IMBH. 

Taking into account the unknown location of the black hole, we examined a field of 
~ ±0.5" in right ascension and dechnation around the center of IRS 13E with a grid step of 
0.01". The residual proper motion of each star results in a required lower mass for the IMBH 
- dependent on its location - in order to bind the star to IRS 13E. Thus, one can obtain 
four maps for the sources El, E2, E3, and E4. The four maps were combined by taking 
the maximum (otherwise the system would not be bound) of the calculated mass limits at a 
given location. We show the resulting map of lower mass limits in the right panel of Fig. 2. 
The significance of the calculated masses over the field is > 4 a. Since El has the proper 
motion with the strongest deviation from the mean proper motion of the ensemble, it is the 
star that gives the strongest constraints. There is a narrow minimum in the mass map for 
black hole locations close to El (~0.1" NW of El). In this most conservative case, the mass 
of the black hole would have to be larger than 7000 M© (with a significance of 4 a). If the 
black hole were located between E2 and E4, the location used by Maillard et al. (2004) for 
their mass estimate, it would have to have a mass > 15000 M© (significance > lOcr). 

These are only lower limits that rely on projected velocities and conservative estimates. 
If IRS 13E is indeed a bound system, a more realistic mass might be derived with the LM 
estimator. Again, there is a minimum near El for the black hole mass. However, even in 
this conservative case, the LM estimator gives a mass of ~50000±15000 M©. 

If we include the source E6 in the above described analysis, the required black hole 
masses would at least double, due to the increased velocity gradient across the complex. 
Hence, we conclude that it is safe to assume that the required mass to bind the members 
of the IRS 13E complex must be at least 10"^ Mq. Assuming the extreme case that the 4 
brightest stars in IRS 13E are similarly massive as the binary IRS 16SW (>100Mq, sec Ott 
et al. 1999), it appears conservative to assume that the stellar mass in IRS 13E does not 
exceed 10^ Mq and can thus be neglected in this analysis. 

4. Discussion 

IRS 13E is located at a projected distance of ~3.5" or ~130mpc from the 3.5 x 
10^ Mq black hole Sgr A*. It may be associated with the counter-rotating disk of young 
stars (Genzel et al. 2003). In this case, it should be located at ~120mpc behind the plane 
of the sky. With a radius of ~0.25", a mass of roughly 1000 Mq then would be sufficient 
to protect the system from tidal disruption. However, as the analysis above shows, the real 
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constraints for binding IRS 13E gravitationally result from the intrinsic proper motions of 
the sources in this complex. We took into account the unknown location of the hypothetical 
black hole and various possibilities for the systemic motion of IRS 13E in the gravitational 
potential of Sgr A*. In the most conservative case, the minimum mass to bind IRS 13E was 
found to be 7000 ± 1800 Mq (with a significance of A a). This would confine the IMBH to 
a narrow region ~ 0.1" NW of El. Outside this region the required mass is greater than 
10^ M0. The LM mass estimator results in even higher masses. 

It is generally accepted that accreting black holes are associated with non-thermal radio 
and X-ray emission. One may argue that an IMBH in IRS 13E could be "starved" , similar 
to Sgr A* (see, e.g., Mclia & Falckc 2001), with its emission below the detection limit of 
current telescopes. However, there are several pieces of evidence that the IRS 13E cluster is 
closely associated with the gas and dust of the mini-spiral: First, there is the close spatial 
relation between IRS 13E and emission from warm dust (see, e.g., the high- resolution AO 
3.8 ;um images in Eckart et al. 2004). Paumard et al. (2004) present an analysis of the 
mini-spiral that also supports a close relation between IRS 13E and the ISM. Also, the stars 
in IRS 13E, with some of them being WR stars, may present favorable sources of gas and 
dust. Moultaka et al. (2004) and Moultaka et al. (2005) present spectroscopic evidence for 
absorption/emission due to gas and dust intrinsic to IRS 13E and for interaction between 
winds from IRS 13E and the surrounding ISM. Hence, the situation for a hypothetical black 
hole in IRS 13E appears to be different from Sgr A*, with plenty of material available for 
accretion. 

As for mm/radio emission, Eckart et al. (2004) show that the 13 mm point source ob- 
served by Zhao & Goss (1998) in the IRS 13E cluster is most likely due to thermal emission 
from the 4//m-excess source E3c in the center of IRS 13E, a dusty WR star (see also Moul- 
taka et al. 2005). The two components E3A and E3B are also proposed as dusty WR stars 
in the analysis by Maillard et al. (2004). As for X-ray emission, Baganoff et al. (2003) re- 
port an X-ray source that they associate with IRS 13E. Using the position of this source 
and of Sgr A* as listed in their Table 1, we located this X-ray source in the NIR reference 
frame. Its position is marked with a 1 o" error box in Fig. 1. The X-ray source appears offset 
from IRS 13E as also discussed by Maillard et al. (2004) . Coker et al. (2002) interpret the 
spectrum of the source as consistent with that of a post-LBV WR colliding wind binary 
(although none of the bright IRS 13E members appears to be associated with the source, 
see Fig. 1). As for a very recent analysis of this X-ray source, it appears to be located a few 
0.1" further north than indicated in Fig. 1, so that E5 hes within the positional uncertainty. 
No variability of the source was found and it appears to be explained well by colliding winds 
(F. Baganoff, priv. comm.). Therefore there appears to be no compelling evidence for an 
accreting IMBH in IRS 13E. Furthermore, spectroscopy at 4 /im (Moultaka et al. 2004) and 
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5 //m (Moultaka et al. 2005) does not reveal any strong broad lines or emission from highly 
excited ions as they might be expected near an accreting black hole. 

Infall and dissolution of a massive cluster into the central parsec seems an attractive 
hypothesis for delivering young stars to the region near Sgr A*. Generally, the conditions 
required for a cluster to form near the GC and to fall in towards the central parsec to deposit 
young, massive stars there within their lifetime, are extreme, i.e., high core densities and 
large cluster masses (see discussions in Gerhard 2001; Genzel et al. 2003; Kim & Morris 2003; 
Kim, Figer, & Morris 2004; Maillard et al. 2004). An IMBH in the core of such a cluster 
was suggested by Hansen & Milosavljevic (2003) as a means of increasing the efficiency of 
dynamical friction on the cluster, and thus to relax somewhat the demands on its total mass 
and core density. Kim, Figer, & Morris (2004), however, conclude that an IMBH in the 
cluster core must comprise an unrealistic ~10% of the total cluster mass in order for this 
mechanism to be effective. 

From the proper motions presented here, the required minimum mass for a hypothetical 
IMBH in order to bind the members of IRS 13E was derived to be > 10"^ Mq. Portegies-Zwart 

6 McMillan (2002) analyze the growth of an IMBH by runaway coUisions in a dense stellar 
cluster and conclude that such an object may form in dense clusters and contain of the order 
0.1% of the total stellar mass. Given the high mass estimate for an IMBH in IRS 13E derived 
in this work, this would mean an unrealistically high mass of > 10^ Mq for the progenitor 
cluster. In this context it is also important to point out the results of Reid & Brunthaler 
(2004) that exclude the existence of a secondary black hole with masses > W^Mq and semi- 
major axes between 10^ and 10^ AU (corresponding to ~0.1" to 12.5" angular distance from 
Sgr A*). We cannot yet exclude an IMBH of up to a few thousand solar masses in IRS 13E, 
but our results suggest that in this case IRS 13E must be in the process of dissolution. In 
this context, we would like to point out the recent work by Levin, Wu and Thommes (astro- 
ph/0502143). They have simulated the infall and dissolution of a stellar cluster with an 
IMBH in the GC and concluded that a hypothetical IMBH that binds the IRS 13E complex 
could not have delivered all of the young stars in the GG to their present location. 

The positions and proper motions of the IRS 13E stars are consistent with their being 
part of the counter-clockwise disk/ring of young stars (Genzel et al. 2003, see also T. Paumard 
et al., 2005, in preparation). As for the unusual clustering of the stars in IRS 13E, the 
complex may be, on the one hand, a (quasi-)permanent feature, i.e. gravitationally bound 
by a hypothetical IMBH. Our analysis shows that in this case, if the star El is indeed part 
of the complex, the IMBH would have to be unsually massive. This makes its existence 
questionable because of the required large mass of the progenitor cluster and because of the 
lack of a clear identification of a non-thermal, variable X-ray source centered within IRS 13E. 
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However, we cannot completely discard the IMBH hypothesis. On the other hand, IRS 13E 
may be a temporary feature, i.e., either a chance association (which is appears somewhat 
unlikely, even if the stars are part of a disk/ring) or a cluster in the process of dissolution. 
At this point, we would like to point out that an unusual concentration of co-moving stars 
has recently also been identified among the members of the IRS 16 complex in the clockwise 
rotating disk of stars (see Lu 2004). Detailed future studies of the proper motions, line-of- 
sight velocities, distribution, and spectral types of the stars in the central parsec of the GC 
nuclear cluster are clearly needed in order to understand better the phase-space clustering 
of the stars in the IRS 13E and IRS 16 complexes. 

Based on observations obtained at the Gemini Observatory, which is operated by the 
Association of Universities for Research in Astronomy, Inc., under a cooperative agree- 
ment with the NSF on behalf of the Gemini partnership: the National Science Foundation 

(United States), the Particle Physics and Astronomy Research Council (United Kingdom), 
the National Research Council ( Canada), CONICYT ( Chile), the Australian Research Coun- 
cil (Australia), CNPq (Brazil) and CONICET (Argentina). 
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Offset from SgrA*, RA 

Fig. 1. — Proper motions of stars in the vicinity of IRS 13E. The underlying image is a LR- 
deconvolved, beam-restored (60 mas FWHM Gaussian beam) NACO K-band image from 8 
July 2004. Two arrows are shown for each source, indicating the ±3 a uncertainty of the 
direction of its proper motion. The lengths of the arrows correspond to the magnitude of 
the respective velocity. The dashed box marked with an "X" designates the position (la 
uncertainty) of the X-ray source near IRS 13E reported by Baganoff et al. (2003). 
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Fig. 2. — Left panel: Zoom into the image shown in Fig. 1. Shown are the residual proper 
motions of the IRS 13E sources after subtraction of the average projected velocity (indicated 
by the dashed arrows) of the sources El, E2, E3, and E4. The opening angles define the la 
uncertainty of the directions. The superposed light grey arrows give the magnitude of the 
velocities after subtraction of a 1 a error. Right panel: Minimum mass required for an IMBH 
in order to gravitationally bind IRS 13E. This mass depends on the location of the IMBH 
and was calculated with the residual velocities of the IRS 13E sources. Contours at 7, 8, 9, 
10, 15, 20, and SOxIO^Mq, growing larger toward lighter shades of gray. The significance of 
the masses is 4-16 a. 
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Table 1: Positions and velocities of sources in IRS 13E (labeling see Fig. 1). Positions are 
given as offsets from Sgr A* in arc-seconds at the epoch 2004.73. 

Source ID RA ["] Dec ["] v^a [kms"!] v^ec [kms"!] 



El 


-2.961 


-1.645 


-206±4 


-51±4 


E2 


-3.171 


-1.732 


-311±4 


77±8 


E3 


-3.191 


-1.527 


-265±17 


171±11 


E4 


-3.210 


-1.411 


-210±9 


124±16 


E6 


-3.013 


-1.422 


-128±16 


41 ±13 



